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The primordial spectrum of cosmological tensor perturbations is considered as a possible probe 
of quantum gravity effects. Together with string theory, loop quantum gravity is one of the most 
promising frameworks to study quantum effects in the early universe. We show that the associated 
corrections should modify the potential seen by gravitational waves during the inflationary amplifi- 
^ , 1 cation. The resulting power spectrum should exhibit a characteristic tilt. This opens a new window 

for cosmological tests of quantum gravity. 
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\l ■ Introduction- A fully convincing quantum theory of gravity is still missing. Beyond the internal theoretical 
difficulties, experimental probes are extremely difficult to find. The two main paradigms to investigate quantum 
gravitational effects are unquestionably string theory and Loop Quantum Gravity (LQG). This latter approach is a 
' promising way to derive a background independent quantum theory of spacetime based only on the experimentally 
. well confirmed principles of general relativity (GR) and quantum mechanics. Introductions can be found in Rovelli's 
book [l| and Smolin's review article 0. Based on the reformulation of GR as a kind of gauge theory obtained 
by Sen 0] and Ashtekar 0] in the 1980's, LQG is now a language and a dynamical framework which leads to a 
mathematically coherent description of the physics of quantum spacetime. Important problems are still to be solved 
in LQG, starting with an explicit semi-classical limit recovering GR, but the important results obtained both in black 
hole physics (see, e.g., and cosmology (see, e.g., Q) urge an experimental test of the theory. There are predictions 
\ for observable Planck scale deviations from energy momentum relations Q that lead to possibly measurable effects 
for and very high energy gamma-ray experiments. Nothing conclusive has however yet emerged from the associated 
analysis. In this article, we focus on a new approach to search for observational probes of LQG [1]. ft consists 
in investigating into the details the way the amplification occurs during the inflationary phase when taking into 
account the modified dispersion relation describing the propagation of gravitons. First, the general framework is 
described together with the Shrodinger-like equation of motion which encodes the LQG corrections in the potential 
term. Then, the primordial tensor power spectrum is explicitly derived and the tilt is given as a function of the 
LQG Barbero-fmmirzi parameter for a specific and favoured value of the n parameter. Both de Sitter and power-law 
inflations are considered. Finally, the validity of the overall scheme is checked and a few points are made about 
possible developments. 

Two basic types of quantum corrections are expected from the Hamiltonian of LQG. These corrections arise from 
inverse powers of the densitized triad and from the fact that loop quantization is based on holonomies, i.e. exponentials 



OS 
O 



X 



•Electronic address: 
t Electronic address: 



grain@apc.univ-paris7.ft'] 
aurelien.barrau@cern.chr 



2 



of the connection, rather than direct connection components. In the following article we focus on holonomy corrections 
as the inverse volume terms were found to exhibit a fiducial cell dependence and not to be adequate to describe models 
with a flat background Q. Unlike studies based on the controversial superinflation scenario (see, e.g., we choose 
in this article to focus on a conservative approach where the background dynamics is supposed to be described by 
a standard de-Sitter or power-law inflation whereas the LQG terms act as corrections on the mode propagation. 
This picture is meaningful for at least two reasons. First, at the fundamental level, because it is highly probable 
PH (although still questioned that superinflation cannot solve alone the problems of cosmology and should be 
considered as an appropriate way to set the initial conditions for the standard inflation. Our approach holds in this 
second phase. Then, at the heuristic level, because in such an intricate framework, it is very useful not to mix all the 
effects. A full LQG treatment of the cosm oloe ical evolution is not yet possible and the aim of this paper is to focus 
on a well determined effect. It was shown |13j that the propagation of gravitons in a FLRW universe is given, when 
holonomy corrections (which provide higher order and higher spatial derivative terms) are taken into account, by the 
following equation of motion: 
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where rj is the conformal time defined by drj ~ dt/a{t), p is a parameter related to the action of the fundamental 
Hamiltonian on a lattice state that can be understood as the coordinate size of a loop whose holonomy is used to 
quantize the Ashtekar curvature components, n is so that p depends on the triad component through p ^ p"', and 
7 is the Barbero-Immirzi parameter. The right-hand side of this differential equation corresponds to the source 
term of gravitational radiations. It also receives corrections from holonomies and vanishes in the absence of matter. 
The friction term and the last term of the left-hand side are given by the background evolution, as solved in the 
LQG framework. One can then compute the equation of motion for tensor perturbation modes with fi — (p/A) , 



n G [—1/2,0]. The value of n depends on the scheme adopted to quantize holonomies. Furthermore |13| . p 



and A has to be chosen so that p has the dimension of a length. The exact value of A and its dependence upon n are 
still under debate. For the specific case n = —1/2, it was shown Q that A = 2^/3Trj£^i and it seems quite natural to 
phenomenologically parametrize phy p = aipip". Using Eq. (29) to (31) from Bojowald & Hossain [ij], the equation 
of motion ([T]) can be re-written as a function of the cosmological parameters (the scale factor 0(77) and the energy 
density of the background p{r])) and of three LQG parameters (n, a and 7): 
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where we have replace ip\ by 1/A/pi and 5* ~ IGnGSl. Introducing a new field $^ — a{rj)h\, Eq. ([2]) now reads: 
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where d is the second derivative according to the conformal time 77. Decomposing the field into its spatial Fourier 
modes, $(a;) = J d^k{27r)''^(t)k exp {ik ■ x), and using the standard Fricdmann equation relating the Hubble constant 
to the energy density of the background, one has to deal with the following equation: 
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It should be underlined that strictly speaking the tree-level Friedman equation in LQG reads as 
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where pcrit ^ ppi = ^/Gl%i in 01 ■ Our whole approach consistently assumes p <^ Pcrit- 

Eq. (|4]) is easily interpreted if one remembers the dynamical equation for gravitons in a FLRW background without 
LQG corrections: 
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with SI the source term in GR. Holonomy corrections appear as a modification of the dispersion relation. This 
modification is encoded in the last term of the left-hand side of Eq. ^ and depends on the dynamics of the 
universe through the scale factor, on its content through the energy density of the background and on the LQG 
parameters. It can also be noticed that, in addition to its time dependence through the scale factor, the correc- 
tion term scales, as expected, as i?background/-Mpi- This equation can be interpreted as a Schrodinger equation, 
being the energy and d/a + 2naj'^ (ct/a^)* a'*+'*" being the potential V. In GR, the potential is simply given by d/a. 



Power spectra- To derive explicitly the power spectrum, we focus on the case n = —1/2. There is a lively debate 
within the LQG community between two versions of the theory and we begin by considering the Ashtekar, Pawlowski 
and Singh scheme. In this approach, the choice n = —1/2 is necessary (see, e.g., [3]) to ensure that the theory 
is independent on arbitrary background structures introduced in the process of quantization. As we deal with the 
amplification of vacuum quantum fluctuations, the source term in Eq. ^ is set to zero. For a de Sitter (dS) inflation, 
the scale factor is given by 0(77) = —l/Hrj with 77 < 0. Therefore, the equation of motion for gravitons 

' ' k U>fc = ; v - - = 2- q;' 



can be solved by a Linear Combination (LC) of Hankcl functions [2^ : 

<Pk{l) = V^[AkH,{^k7^) + BkHli-ki^)^ , (6) 

A), and Bk being two constants of integration determined by the initial conditions found by studying the region where 
the adiabatic vacuum can be defined (see, e.g., p^). In this so-called Bunch-Davies vacuum, the solutions are plane 
waves 

lim (^fc(77) = — -= (7) 

kri^-00 A7pi ^2k 

whose matching with the previously given general solution allows to determine Ak and Bk. The choice of the vacuum 
is clearly a non-trivial question in this framework. However, to remain consistent with our hypothesis, we assume 
the de-Sitter background to be a correct approximation from the beginning of inflation. In this case, the appropriate 
conditions will inevitably be met in the remote past , i.e. when krj —00. This will not be true anymore when 
considering inverse volume corrections p^ . Finally, the solution is expanded around the high amplification regime, 
which leads to the following power spectrum, defined as PT{k) = 2fc^/(7r^) x \4>k{rjf)\ /a^{rjf) (using the convention 
of (H): 

V^{k)=A^k^-^^ ; At^(^ ' JiJ ) ■ 

It can be seen that LQG corrections do not only modify the normalization of the spectrum but also, and more 
importantly, lead to a departure from scale invariance. The tensor index is given by; 

"T = 3 - 2.. ^ ^ (^) ' + O(i^VMpi) (9) 

which shows that the spectrum becomes blue for a real-valued Barbero-Immirzi parameter (which is the more realistic 
case, as inferred, 6.17., from the entropy of black holes) and red if 7 is imaginary. This can be easily understood at the 
intuitive level. The amplification starts at the critical time rjc = —^/vjk. If 7 e M this is later than in the standard 
general relativistic case and modes are therefore less amplified. As the difference between the value of 77c with LQG 
corrections and without LQG corrections becomes smaller as k increases, short wavelengths will undergo a smaller 
suppression. The spectrum will therefore be blue. 

Several cases can be investigated into more details. First, it is worth considering a power-law inflation. For this 
inflationary model, the scale factor is given by 0(77) — io l??!'^'''^, with /3 < —2. Using S = y/a{f3+ 1)'^jHq/ |/3 + 2| A/pi, 
where Hq = i^^, the potential term now reads 

/3(/3 + l) d^\/3 + 2\' 
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FIG. 1: Primordial power spectrum, in Planck units, for power- law inflation with (3 = -2.5 and Ho/Mpi = 10""*, 10"^ and 
lO^'^ (from the left to the right): dashed lines are for the numerical results, blue lines are for the IR limit and red lines for 
the UV limit of the power spectrum. The UV approximation coincides with the primordial spectrum obtained without LQG 
correction. 



Because P < —2, the LQG corrections dominate at the beginning of inflation and become negligible at the end. 
Moreover, it can be seen from Fig. [5] that adding LQG corrections lowers the effective potential. As the amplification 
occurs when < V{r]), it is obvious from a simple WKB reasoning that high fc-modcs will not be much affected by 
LQG corrections whereas low fc-modes will be less amplified than in the standard power-law case. LQG corrections 
being subdominant at the end of inflation, the mode functions are still well approximated by a LC of Hankel functions 
like in Eq. ^ with v = \l3 + 1/2|. The primordial power spectrum is then given as a function of Ak and Bk by taking 
the asymptotic expansion —krj—^Ooi the Hankel functions: 

Vrik) = H'o (^^^^^^^^ k'-'-'+'lAk-B.f. (11) 

With the considered potential, the equation of motion for gravitons cannot be solved analytically, though the asymp- 
totic limits can be found. Because the LQG term is subdominant at the end of inflation, the ampliflcation is effective 
only in the UV regime (i.e. k —f oo) and the power spectrum should be approximated by the standard prediction of 
power-law inflation for high values of the wavenumber. This leads to 

Then, it is possible to compute the spectrum in the IR limit (i.e. k —f 0). The fc— dependant solution given in Eq. 

valid at the end of inflation, is matched to a solution valid during the entire inflationnary era, but only for fc = 0. 
By setting fc = and with the change of variables —77 = exp {x/ |/3 + 2|) , 0o = u{x) exp{x/2 \P + 2|), the equation 
of motion takes the form of a Bessel equation. This allows us to find a solution for the (fc = 0)-modes function, 
valid throughout the whole inflationary era, as a LC of Hankel functions. The coefficients of this solution are found 
requiring the Wronskian to be equal to W ^ —ildn/Mpi j23j : 

Mv) - ^^VW+^\H, {S\vt'^'^) , (13) 

with fi — \P + l/2\ / \P + 2\. Finally, the matching between the two solutions is performed at the end of inflation by 
taking the asymptotic limit of the solution of ([6]) when — fc?; and the one of solution (fT3|l when — ?y 0. This 
leads to: 

\A. -B,f^ ' 1/3 + 21 fc-^-. (14) 
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This allows us to derive the primordial power spectrum in the IR and UV regimes: 

In both regimes, the power spectrum is well described by a power-law with a red (blue) spectrum in the UV (IR) 
region. To compute the power spectrum throughout the entire range of wavenumbers, the differential equation has 
been numerically solved with the Bunch-Davies vacuum as initial conditions. The results are displayed on Fig. [1] 
This shows the very good agreement between the numerical computation and the UV and IR approximations. The 
value of (5 was deliberately chosen very small to make the reading easier. As the UV approximation is also the 
primordial spectrum obtained without LQG corrections, it can be seen that adding the LQG correction term leads to 
a suppression of the long wavelength perturbations as compared to the prediction of GR. 

The results obtained in the power-law inflation framework can be directly applied to slow-roll inflation, which is a 
realistic scenario in agreement with WMAP data. In particular, the asymptotic expansions of Eqs. (15) and (16) are 
still valid with the appropriate replacement of/3-|-lby— 1 — e [23|. This leads to: 



^t(A:)=(^^) [l-e(2-ln2 + 21n(-fc77,))]r2 Q-fe) fc|Afc-Bfc^ (17) 

where iJg has been expressed as a function of rjc and iJc, respectively the conformal time and Hubble constant at 
horizon crossing . A Taylor expansion in e has also been performed. 

Prospects and validity of the approach- In the Ashtekar, Pawlowski and Singh approach considered so far, n must 
be equal to -1/2. However, Bojowald has criticized this conclusion (see, e.g., [l^l for a very recent discussion) and, for 
the completeness of our conclusions, it worth investigating the case n > —1/2 at the qualitative level. For the three 
inflationary models studied in this article (dS, power-law and slow-roll), the potential reads as V{rj) = Ai/77^ — A2/ \rif, 
where A, are positive constants such that A2 ^ Ai and p = 4 — 4n(/3 4- 1). Depending on the values of n and /?, the 
exponent of the LQG correction terms p is either smaller or greater than 2. (When n = —1/2, p is always smaller 
than 2). The shape of the potential is displayed on Fig. [2] together with the GR potential. The case p < 2 (dashed 
line) was studied in the previous section, where numerical investigations confirmed the naive WKB-based idea that 
amplification is less effective in the IR regime while roughly the same than in GR in the UV regime. The situation 
drastically differs when p > 2 (dotted line) as the LQG term now dominates at the end of inflation (this leads to 
the drop of the potential when — ?/ 0). With the same WKB reasoning, it is clear that the amplification is now 
less effective for high values of k. Moreover, because the potential admits a maximum, the amplification totally 
disappears for k 00. Finally, as the potential becomes negative valued, the mode functions can start to re-oscillate 
at the very end of inflation. This fundamentally changes the standard picture of the generation of perturbations as 
the temporal coherence of the two-point correlation function is restored before the transition from inflation to the 
radiation-dominated era. Even with (3 very close to -2, as expected from WMAP data [l3|, this can happen for n 
slightly above -1/2. 

It is well known that a complete description of the primordial fluctuations amplified during inflation requires, in 
principle, to derive the full quantum statistical properties of the Bunch-Davies vacuum (the in vacuum) in terms 
of quanta defined in the radiation-dominated era (the out quanta). In the standard GR case, the quantum statis- 
tical distribution is maximally coherent due to the high enough amplification of quantum fluctuations (i.e. when 
k/aH 0). This high level of coherence ensures the classical description to be a good approximation at the level 
of power spectra. In other words tensor fluctuations can be described by a classical, Gaussian distribution with zero 
mean and a variance proportional to the power spectrum (see [l^ for a detailed discussion of the link between the 
apparent classicality and the coherence of the quantum distribution). When LQG corrections are considered, the 
above statement still holds when n = — 1/2 simply because LQG corrections are negligible at the end of inflation. 
When n > —1/2, LQG correction are now dominant at the end of inflation and a complete analysis is needed. For 
super-horizon modes, an asymptotic solution in terms of Bessel functions of the fleld equation can be found, allowing 
us to explicitly compute the Bogoliubov coefficients relating the in vacuum to the out quanta. In the limit k/aH — > 0, 
the two Bogoliubov coefficients are related one to the other by 



Pk = ai exp {--2ik/aH). 



(18) 
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FIG. 2: Qualitative shape of the effective potential during inflation without LQG corrections (solid line) and with LQG 
corrections: the dashed line is for p < 2 and the dotted one for p > 2. 



Because k/aH <^ 1, the above-mentioned relation ensures that the distribution is maximally coherent and the 
classical description is still a good approximation. In other vifords, the classical description valid in a general 
relativistic framework still holds when LQG settings are considered i/ inflation lasts long enough. 

Conclusion- Testing quantum theories of gravity is probably one of the most important challenge of current fun- 
damental physics. Loop Quantum Gravity has not yet been experimentally probed but it seems that cosmological 
observations could allow for a clear signature of LQG effects. This article opens a new window on quantum gravity 
by suggesting a way to possibly observe holonomy corrections. Although B-mode detection could be achieved by 
the Planck satellite and is the main goal of several dedicated experiments for the forthcoming decade, the signal 
remains very difficult to extract from the lensing background and some further refinements are required to quantify 
the amplitude of the expected LQG effects. Furthermore, string gaz cosmology also predicts a deviation from scale- 
invariance [20| and the discriminating criteria should be found. Most importantly, the formalism established in this 
article should be used for LQG corrections to scalar perturbations that are just being investigated [2l| and could be 
even more promising from the observational viewpoint. 
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